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We calculate the upper critical field in superconductors without inversion symmetry at arbitrary 
temperatures, in the presence of scalar impurities. Both orbital and spin (paramagnetic) mechanisms 
of pair breaking are considered. The superconducting phase transition at nonzero field occurs into 
a helical vortex state, in which the order parameter is modulated along the applied field (in a cubic 
crystal). The helical state is stable with respect to disorder. However, if the difference between 
the densities of states in the electron bands split by spin-orbit coupling is neglected, then the order 
parameter modulation is present only at low temperatures, and only if the system is sufficiently clean 
and paramagnetically limited. This state resembles the Larkin-Ovchinnikov-Fulde-Ferrell state in 
centrosymmetric superconductors. 



PACS numbers: 74.20.-z, 74.25. Op, 74.25.Dw 



I. INTRODUCTION 

Recently, superconductivity has been discovered in 
a number of noncentrosymmetric compounds. Physi- 
cal properties of these materials vary considerably, from 
CePtaSi (Ref. H, in which strong electron correlations 
are responsible for a heavy-fermion behavior, and the 
superconductivity is likely anisotropic with gap nodes, 
to LizPdaB and LisPtsB (Ref. [1), which show rather 
conventional features, described by the Bardeen-Cooper- 
Schrieffer (BCS) theory of phonon-mediated pairing. 

The absence of inversion symmetry in the crys- 
tal lattice brings about important qualitative changes, 
both in normal and superconducting properties, com- 
pared with the centrosymmetric case. These dif- 
ferences are highlighted, in particular, by different 
responses to an external magnetic field. Distinc- 
tive features of noncentrosymmetric superconductors in- 
clude a strongly anisotropic spin susceptibility with 
a large residual component ^i^i^i^iiiSii^ magnetoelectric 
effect ^iiSiiiiii^ and novel nonuniform ("helical") super- 
conducting statesJ^ii^ii^ 

In this paper, we extend the classic theory of the up- 
per critical field, Hc2{T), in BCS superconductors, which 
was developed in Refs. 16 and [itI, see also Ref. 18, to 
the noncentrosymmetric case. The spin-orbit (SO) cou- 
pling of electrons with a noncentrosymmetric crystal lat- 
tice considerably changes the nature of single-electron 
states, lifting spin degeneracy of the energy bands. This 
modifies the Zeeman coupling of band electrons with 
magnetic field, which in turn affects the paramagnetic 
pair breaking. The magnetic phase diagram of noncen- 
trosymmetric superconductors has been previously stud- 
ied in a two-dimensional case ;^^'^" in which the orbital 
effects are not important. In Ref. the upper crit- 
ical field for a three-dimensional Rashba superconduc- 
tor was calculated in the absence of impurities. The ef- 
fects of disorder on Hc2 in the Ginzburg-Landau regime 
near the zero-field critical temperature have been inves- 
tigated in Ref. [2l| for an arbitrary pairing symmetry, but 
without paramagnetic effects. In this paper, we include 



both orbital and paramagnetic mechanisms of pair break- 
ing, as well as scalar disorder, at arbitrary temperatures. 
We consider a "minimal" model of a three-dimensional 
weak-coupling superconductor of cubic symmetry, with 
an isotropic parabolic band split in two by the SO cou- 
pling. We assume that the ratio of the SO band splitting 
to the Fermi energy is small (which is the case in most 
known noncentrosymmetric materials), and neglect the 
impurity-induced mixing of the singlet and triplet com- 
ponents of the gap function. 

The paper is organized as follows: In Sec. [Til we derive 
the equations for Hc2 (T) , relegating some of the technical 
details to Appendix. In Sec. IIIIl the transition tempera- 
ture as a function of the external field is calculated in the 
Ginzburg-Landau region near the zero-field critical tem- 
perature TcQ. In Sec. llVi we consider the purely para- 
magnetic limit, both in the clean and disordered cases. 
The general case, with the orbital pair breaking included, 
is studied in Sec. |Vl Sec. IVII contains a discussion of our 
results. Throughout the paper we use the units in which 
h = ks — ^■ 



II. BASIC EQUATIONS 

Let us consider a noncentrosymmetric superconductor 
with the Hamiltonian given hy H = Hq + Himp + -ffmt ■ 
The first term, 



(1) 



describes non-interacting electrons in the crystal lattice 
potential, where a,/3 =t,i are spin indices, eo(fc) is the 
quasiparticle energy, and & are the Pauli matrices. We 
assume a parabolic band and include the chemical po- 
tential in the dispersion function: eo{k) = k'^ /2m* — ep, 
where m* is the effective mass, ep — kQ/2m*, and fcg 
is the Fermi wave vector in the absence of the SO cou- 
pling. In Eq. IT]) and everywhere below, summation over 
repeated spin indices is implied, while summation over 
band indices is always shown explicitly. 
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The second term in Eq. ([T]) describes a Rashba-type 
SO coupling of electrons with the crystal lattice^^^ We 
focus on the case of a noncentrosymmetric cubic crys- 
tal with the point group Gr — O, which is applicable, 
for instance, to Li2(Pdi_2;,Pta;)3B. The simplest expres- 
sion for the SO coupling compatible with all symmetry 
requirements has the following form: 



7(fc) = 7ofc, 



(2) 



where 70 is a constant. It is convenient to characterize 
the SO coupling strength by a dimensionless parameter 



Q = 



77i*|7o| 

kn 



(3) 



Diagonalization of the Hamiltonian ^ yields two non- 
degenerate bands labelled by helicity A = ±, which are 
described by the following dispersion functions: 



CA(fe) = eo(fe) + A|7(fc)| 



2m 



M + A|7o|fc. (4) 



The SO band splitting is given by Eso — 2|7o|fco, there- 
fore g = Eso/4:eF- While the two Fermi surfaces, de- 
fined by the equations ^\{k) = 0, have different radii: 
kp.x — fco ( V 1 + ^ ^q)} the Fermi velocities are the 
same: v\ — vpk, where vp — {ko/m*)\/l + g^. 

Scattering of electrons at isotropic scalar impurities is 
introduced according to 



Hi, 



SrU{r)i^i{r)^^{r), 



(5) 



where the impurity potential U (r) is a random func- 
tion with zero mean and the correlator {U{ri)U{r2)) = 
nimpU^5(ri — r2), Ji-imp is the impurity concentration, 
and Uq has the meaning of the strength of an individual 
point-like impurity. The field operators have the usual 
form: '0q('") — V~^/^ e^^^aka, where V is the system 
volume. 

We assume that there is a local attraction between 
electrons, e.g. due to phonons, and describe the Cooper 
pairing by a BCS-likc Hamiltonian: 

H,^t =-V f d3r^|(r)7A|(r)Vi(r)^t(r), (6) 



where 1^ > is the coupling constant. A detailed analysis 
of the relation between the microscopic pairing interac- 
tion and the gap symmetry in the band representation 
can be found in Ref. [H. In particular, in the BCS- 
contact model (|6]) the local pairing of electrons with op- 
posite spins translates into the same-helicity pairing, the 
superconducting gap function has only intraband compo- 
nents, and the order parameter is represented by a single 
complex function ri{r). 

External magnetic field can be included in the elec- 
tron band theory by making the so-called Peierls 
substitution^^ in the dispersion functions (,\{k): 



K 



.d_ 
dr 



-Mr). 



(7) 



where e is the absolute value of the electron charge 
and A(r) is the magnetic vector potential. Near the 
upper critical field, the magnetic induction is uniform, 
B{r) — H, so that A{r) — {H x r)/2 in the symmet- 
ric gauge. This approach has been used to microscopi- 
cally derive the Ginzburg-Landau functional for noncen- 
trosymmetric superconductors, both in the clean^^ and 
disordered^- cases. However, the bands (jH) are nonana- 
lytic in k, and the Peierls substitution leads to ill-defined 
operators. To avoid this, one can make the substitution 
([7]) directly in the original Hamiltonian [Eq. Then, 
the SO coupling, the impurity scattering, and the mag- 
netic field are all incorporated in the following single- 
electron Hamiltonian in the coordinate-spin representa- 
tion: 

q 

h^ — +-faK& + ^fiBH& + U{r)^eF, (8) 

where g is the Lande factor, and /is is the Bohr magne- 
ton. The final expressions for observable properties, in 
particular the upper critical field, do not actually depend 
on whether the Peierls substitution is made in the spin 
or band representations. 



A. Gap equations 

Following the standard textbook procedure, one can 
calculate the difference between the disorder-averaged 
free energies in the superconducting and normal states, 
at the same temperature and field. In the vicinity of the 
upper critical field Hc2{T), one can keep only the terms 
quadratic in the order parameter in the free-energy ex- 
pansion: 



^ = j J dridr27f{ri)S{ri,r2)r]{r2) 



(9) 



We assume that the superconducting phase transition is 
of second order. It should be noted that the validity 
of this assumption at all temperatures is not obvious, 
see, e.g., a discussion in Sec. IIV CI below, and should be 
verified by examining the higher-order terms in the free 
energy, which is beyond the scope of the present work. 
The kernel in Eq. ^ has the following form: 

S{r,,r2) = ^d{r, - rs) - 'A(ri, ra; a;„), (10) 



where 



A(ri,r2;w„) = ^gl^g^s 



>i{Gi3j{ri,r2;uJn)Gas{ri,r2; -w,,))^^^, (11) 

ujn = (2n + l)nT is the fermionic Matsubara frequency, 
and g ~ ia2. The prime in the second term in Eq. (jlOp 
means that the summation is limited to |u;„| < Wc, where 
LUc is the BCS frequency cutoff. The angular brackets 
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denote the impurity averaging, and G'(ri,r2;a;„) is the 
Matsubara Green's function of electrons in the normal 
state, which satisfies the equation 



{iuj„ - /ii)G(ri,r2;w„) = d{ri - r2), 



(12) 



with the Hamiltonian h given by expression ([8]). The 
subscript "1" means that h acts on the first argument of 
the Green's function. The critical temperature at a given 
field, or inversely the critical field at a given temperature, 
is found from the condition that the lowest eigenvalue of 
the operator S, defined by the kernel (jlOp. is zero. 

At zero field, Eq. (fT2)l yields the following expression 
for the average Green's function:^^ 

G{k, c^„) = nA(fc)GA(fe, c^„), (13) 



where 



Ilxik) = 



1 + Xk& 



(14) 



are the band projection operators (fc = fc/|fc|), and 



Gx{k,UJn) 



lUJn 



£,x{k) +irsignw„' 



(15) 



are the electron Green's functions in the band represen- 
tation. Here F = 1/2t is the elastic scattering rate, 
T = {2TrnirnpUQ Nf)^^ is the electron mean free time 
due to impurities, Np = (7V_|_ + N_)/2, and Nx is the 
Fermi-level density of states in the Ath band. 

Using the standard techniques,— the impurity average 
of the product of the two Green's functions in Eq. (fTTj) 
can be represented graphically by ladder diagrams, see 
Fig. [T] (we assume the disorder to be sufficiently weak 
for the diagrams with crossed impurity lines to be negli- 
gible). Summation of the diagrams is facilitated by rep- 
resenting the impurity line (which corresponds to each 
"rung" of the ladder) as a sum of spin-singlet and spin- 
triplet terms: 

(16) 

where g — ia&2- We introduce an impurity-renormalized 
gap function D^jsiq, uJn), which satisfies an integral equa- 
tion 



D{q,uJn) = viQ)9 

1 f (Pk 

+ ^n,^pU^g J j^--:^tT[g''G{k + q,ujn) 

xI)(q,w„)G^(-fc,-tj„)] 

1 f (Pk 

+ ^n„npUag J tr [g^G{k + g, w„) 

xD{q, w„)G^(-fc, -LJn)] ■ 



(17) 



g+ g + g+ 



— ^ — 

P o 



g + 



FIG. 1: Impurity ladder diagrams in the Cooper channel. 
Lines with arrows correspond to the average Green's func- 
tions of electrons in the spin representation, g — i(j2, and the 
impurity (dashed) lines are defined in the text, see Eq. (I16p . 



Seeking solution in the form 

b{q, ujn) = do(q, ^n)9 + d{q, Lj„)g, (18) 

we obtain the following equations for da{v^LOn) (a = 
0,1,2,3): 

3 

da{q,uJn) = ■n{q)5aQ + TYyab(,q,^n)db{q,uJn), (19) 

6=0 

where ri{q) is the Fourier-transform of the order param- 
eter, and 



3^afc(g,t^n) 

d^k 



1 



2'kNf 



tr [gtG(fc + q, c.„)g,G^(-fe, -co„)] , (20) 

with go = g, and gi = cji for i = 1, 2, 3. We see that, in 
addition to the spin-singlet component dQ(q,u)n) of the 
gap function, impurity scattering can induce a nonzero 
spin-triplet component d(q,LJn)- Note that the free en- 
ergy depends only on the singlet component: From Eqs. 
(HH) and (ini), we obtain: 



d^q 
(2^ 



v*iq) 



V 



viq) 



-ttNfTY^ 



,do{q,uJn) - 77(9) 



(21) 



Substituting the Green's functions (|T3|) into Eqs. pOI 
and calculating the spin traces, we obtain for the singlet 
singlet part of the 4x4 matrix y-. 



ym{q,^n) 



1 



2'kNf 
d^k 



/d k 
-J^^G\{k + q,ujn)Gx{-k, ~uj„) 



1 

l^nl + r + ifi signtJn / f, 

1 



where 



\u!n \ + T + signwn / ' 
N± 



(22) 



(23) 
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are the fractional densities of states in the two bands, 
V,{k,q) = VFkq/2, and denotes the averaging 

over the directions of k. The parameter 



6 = 



N+ - 



(24) 



which characterizes the difference between the band den- 
sities of states, can be expressed in terms of the SO cou- 
pUng strength ([3]) as follows: 



l + 2g2 ■ 

Similarly, for the singlet-triplet mixing part we obtain 

1 



(25) 



2TTNf 



d?k 
(2^ 



hG\{k + q,LUn)Gx{-k,-uj„) 



|w„| + r + zil signcj„ 



(26) 



Finally, the triplet-triplet part can be represented as fol- 
lows: 



y^J{q,iOn) = yll\q,uj,,)+yl^\q,LO,,), (27) 



where 



1 



2ttNf 



E / j^:^f^ikjGx{k + q,ujn)Gx{~k,~ujn) 



|w„| + r + ifi signa;„ / . ' 

/ k 



and 



(fk 



2ttNf ^ J (27r)^ 



—y 



{^ij kikj iXcijiki) 



^ij kikj iXcijiki 




2^ \ \ujn\+T + + XEso/2)signujn .. ' 

We see that in the band representation the singlet im- 
purity scattering channel, which is described by the first 
term in Eq. (|16p. causes only the scattering of intraband 
pairs between the bands. In contrast, the triplet channel 
can mix intraband and interband pairs, the latter being 

described by y^j . 

The superconducting critical temperature is found by 
setting q — 0, then J^oi = 3^*0 = 0, and the solution of 
Eqs. (fTO|) has the form do — {1 + r/|ti;„|)r;. Substituting 



this into Eq. (|2T|) and comparing the result with Eq. 
we obtain: 



5(, = o)^l-.iv.r^'J- 



(28) 



which yields the superconducting critical temperature: 



^ 2e^ . 
Jco = Wee 

TT 



-1/NfV 



(29) 



where C ~ 0.577 is Euler's constant. Thus there is 
an analog of Anderson's theorem in noncentrosymmet- 
ric superconductors with a BCS-contact pairing interac- 
tion: The zero-field critical temperature is not affected 
by scalar disorder.— 

Now let us turn on a uniform magnetic field H — Hz. 
Its orbital effect is taken into account in the usual way, by 
replacing q ^ D = — iV -I- {2e/c)A in yab{q, w„), which 
become differential operators of infinite order. The gap 
equations then take the following form: 



(1 - Tya^)do - vyoA - ?/, 



(30) 



Solution of these equations in the general case is rather 
cumbersome. In order to make progress, we use the fact 
that in practice the SO coupling is much weaker than the 
Fermi energy. In terms of the parameter (5, this trans- 
lates into the following condition: \5\ <C 1. This allows 
us to neglect the triplet component of the gap function. 
Indeed, using the fact that the singlet-triplet mixing is 
described by J^oi = 3^iOj which are proportional to 5, see 
Eq. (pS)) . one can solve the gap equations ([50]) by iter- 
ations. It is easy to see that d ~ 0{5), therefore the 
correction to do — (1 ~ 1^3^00)^^?? due to the impurity- 
induced singlet-triplet mixing is of the order of 5'^ and 
will be neglected. 



B. Equation for Hc2(T) 

Keeping only the singlet gap component, we have do = 
(1 - rj^oo)"^^- Substituting this in Eq. dH]), we can 
represent the operator S as follows: 

S ^ ^ - ^NfTY^'I^ -^y00{^n)r^%0{^n), (31) 



where 3^oo(wn) is defined by the kernel 



ym{ri,r2;ujn) 



1 



27riVi 



X (G^^(ri, ra; a;„))^^^(G,,5(ri, ra; -c^„))^,^p. (32) 

Although the expression ()31|) is approximate (with the 
corrections of the order of 5^), it has the advantage of 
being relatively simple and captures important physics of 
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the problem, including the properties of various nonuni- 
form superconducting states created by the external field, 
see Sees. IIIHTVl below. The general case of arbitrary SO 
coupling, with both singlet and triplet channels present, 
will be considered elsewhere.-^ 

It is shown in Appendix A that the eigenfunctions of 
3^00 (wn) are the Landau levels \N,p), which are labelled 
by two quantum numbers: a non- negative integer iV and 
a real p. The latter is proportional to the wave vector of 
the order parameter mod ulation along the applied field: 
P — ^Hlzj where in = \/c/eH is the magnetic length. 
The corresponding eigenvalues can be written as follows: 



j)oo(Wn)|A^,P> = yN,pi'^n)\N,p), 



(33) 



where 

Jo 

pi 

X / dsFp{u,s)e-^'^^-''^^^LN[v^{l-s'')]. (34) 
Jo 

Here v — (vp /2£h)u, Ln{x) are the Laguerre polynomi- 
als, 

Fp{u, s) = cos{povs) cos(pws) — S sin(povs) sin(pws),(35) 



Po — g^^BH^H /vp, and 5 is defined by Eq. ([M)) . Note 
that in centra symmetric isotropic superconductors with 
nonmagnetic impurities, the eigenvalues of 3^oo(^^n) are 
given by the same expression p4p . but with Fp{u,s) = 
cos{pov) cos{pvs). This is still different from Eq. (pS)) . 
even if one sets (5 = in the latter. The reason is that 
Eq. (j35p is valid under the assumption that the Zeeman 
energy is small compared with the SO band splitting. 

Since the Landau levels are eigenfunctions of y 00(10,1) 
at all frequencies, the operator S, see Eq. pT|) . is also 
diagonal in the basis of \N,p). Using Eqs. P5|) and (pg)) . 
we can eliminate both the frequency cutoff and the cou- 
pling constant from Eq. (|3ip. In this way we obtain an 
equation for the magnetic field at which a superconduct- 
ing instability characterized by N and p develops: 



In- 



TcQ 



_J yN,pi'^n) 

\u!n\ l-Typf^piuJn) 



(36) 



The upper critical field, Hc2(T), is obtained by maximiz- 
ing the solution of the above equation with respect to 
both A'' and p. 

It is convenient to introduce the reduced temperature, 
magnetic field, and disorder as follows: 

T 2H F 

^=^' ^=7r' C=^, (37) 

where Hq = $o/^Coi '^'o — 7rc/e is the magnetic fiux 
quantum, and ^0 — vp/'^'^Tco is the superconducting co- 
herence length. Eq. ([M]) then takes the form 



i4 = 2E 



n>0 L 



1 



„(JV,Q) 



- t— 

271+1 I 



{t,h) 



(38) 



where 

poo pI 

Jo Jo 



and 



X exp 



^g(p, s) 



L 



N 



cos{ahps) cos{Qps) 

—S sin{ahps) sm{Qps). 



The parameter 



g f^BH( 







2 2'kT, 



cO 



,(39) 



(40) 



(41) 



measures the relative importance of the paramagnetic 
and orbital contributions to the magnetic suppression 
of superconductivity (note that a is proportional to the 
Maki parameter an, introduced in Ref. [27h . For the 
upper critical field in the reduced notations we have 



hc2{t) = max/iAT n(t), 

N,Q 



(42) 



where hN,Q{t) is the solution of Eq. ([38]) . A^ is the Lan- 
dau level index, and Q = ^oQz is the dimensionless wave 
vector of the superconducting instability. 

The purely orbital limit is obtained by formally set- 
ting g to zero. Then a — 0, and Eq. ([55)1 takes 
exactly the same form as in centrosymmetric BCS su- 
perconductors. Therefore, if the Zeeman interaction is 
neglected then the absence of inversion symmetry does 
not bring about any new features in Hc^iT), compared 
with the centrosymmetric case, described by the Helfand- 
Werthamer theory. The maximum critical field corre- 
sponds to A^ = Q = at all temperatures. In particular, 
at T = and in the absence of impurities one obtains: 
Hc2{0) = (e2-C/8)iJo ~ 0.52Ho. This can be used to re- 
late a to experimentally observable quantities as follows: 



0.21 



gc2(0)[r] 



(43) 



Here we assumed g — 2. 



III. GINZBURG-LANDAU REGIME 

At weak external field, the critical temperature is close 
to Tco, i.e. the reduced temperature t is close to 1. In this 
limit, one can solve Eq. psp analytically, by expanding 
^{N,Q) powers of h and Q (we set A^ = 0). We seek 
solution in the form 



tc{h) = l-aih- 02/1^ + 0{h^). 



(44) 



It follows from Eq. (|40|) that the maximum critical tem- 
perature corresponds to 



Q 



-6ah. 



(45) 
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After some straightforward algebra, we obtain the fol- 
fowing expressions for the coefficients in the expansion 
(SI: 



with 



1 2 1 2 1 



(46) 



(47) 



ji>0 



Since our model is valid only in the limit of weak SO 
coupling, \6\ <C 1, we have omitted in these expressions 
the terms containing S^. In the clean limit, we have 



a2 



7C(3) 
12 ' 
49C'(3) 
96 



31C(5) 
40 



7C(3) 2 

a , 

12 



where Ci^) is the Riemann zeta function. 

Since Eq. does not contain a, the slope of the 

upper critical field near Tco is entirely determined by the 
orbital effects. Returning to dimensional variables, we 
obtain: 



dH. 



c2 



dT 



6tT $nr. 



O^cO 



52,1 (0 



(48) 



which has exactly the same form as in isotropic cen- 
trosymmetric superconductors.^'* From this expression it 
follows, in particular, that the upper critical field slope 
in our model is enhanced by disorder, in agreement with 
the results obtained in Ref. Hi]. 

The effects of the Zeeman interaction on the critical 
temperature appear only in the second order in h and are 
described by the last term in Eq. (|47)) . If the difference 
between the band densities of states is taken into account 
(i.e. if (5 ^ 0), then the order parameter is modulated 
along the applied field, with the period given by Eq. ps)) . 

In the limit of large a, the Zeeman term dominates 
and the critical field slope diverges: Hc2{T) ~ \/Tco — T. 
Similar to the orbital critical field ([48]), the paramag- 
netic critical field in the Ginzburg-Landau regime is en- 
hanced by disorder. The paramagnetic pair breaking in 
noncentrosymmetric superconductors has some peculiar 
features, compared with the centrosymmetric case, and 
is studied in detail in the next section. 



IV. PARAMAGNETIC LIMIT 

The purely paramagnetic limit corresponds to a — > cxd 
in Eq. ([55)1 . Due to the fast oscillations of $q, the last 



two factors in the s-integral in Eq. ()39p can be replaced 
by 1, and the solution becomes degenerate with respect 
to N. We have 



n{Q) = \'^P^ 



xRe I 



1 



where 



{2n + l)t + C + i{Q + Xh)k^ 



7 9 fJ-sH 



, (49) 



Calculating the Fermi-surface integrals, we obtain: 
1^- 1 . Q + Xh 



arctan ■ 



2 X + (2n+l)< + C 



(50) 



(51) 



Further analytical progress can be made in two limit- 
ing cases: clean {( = 0) and "dirty" (C ^ 1), while the 
intermediate disorder strengths can only be studied nu- 
merically. 



A. Clean case 

In the clean limit, it is convenient to go back to ex- 
pression ([^5)1 and substitute it into Eq. ([55]) . with the 
following result: 



In- 



11 

1^2 



1 .Q + Xhf 
— h i k 

2 2t 




(52) 



where "^{x) is the digamma function. From this we ob- 
tain: 



t — exp maxX(y, z) 
y 



(53) 



where y = Q/2t, z = h/2t 
'1 



X(y,z) = * 



1^ lra\iiT{l/2 + iy + iXz) 

o z^p^ r~\ ' 

2 V 2/ + 



and T(x) is the Gamma function. At any given z, if the 
maximum of I{y, z) is achieved at y = j/c, then the wave 
vector of the superconducting instability is Q = 2tyi., and 
the corresponding critical field is h — 2tz. The critical 
field of a second-order phase transition into a uniform 
superconducting state can be found by setting y = in 
Eq. ([53]) . In particular, at zero temperature this phase 
transition occurs at ho = l/2e^~"'^ ~ 0.76. 

The pair-breaking effect of the Zeeman interaction can 
be reduced by allowing the pairs to have a nonzero center- 
of-mass momentum. The outcome of the competition 



between the Zeeman energy and the gradient energy de- 
pends on the difference between the densities of states 
in the two bands, which in turn depends on the ratio 
of the SO band sphtting to the Fermi energy, see Eqs. 
([23)1 and ([25]). If it is neglected, i.e. 5 = 0, then one 
obtains from Eq. ([53]) that at z < z* ~ 0.44, which 
corresponds to high temperatures, t > t* ~ 0.68, and 
low fields, h < h* 0.60, the maximum critical field is 
achieved at Q = 0, and the phase transition occurs into 
the uniform superconducting state. However, at i < t* 
the maximum critical field corresponds to Q ^ 0, and the 
phase transition occurs into a nonuniform superconduct- 
ing state, similar to the Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) statci^ The order parameter is modulated along 
the field: r]{r) = ?7ie*^^ -I- ry2e~"^^, where q = Q/Co- The 
coefficients 771^2 are found from the higher-order terms in 
the free energy. In Refs. [13 and [s^, this was done for 
a Rashba superconductor, and it was shown that both 
a multiple-g (or stripe) LOFF state and a single plane 
wave state are possible, depending on temperature. 

If the difference between p+ and p- is taken into ac- 
count, then the right-hand side of Eq. ([5^ is no longer 
even in and the maximum critical field corresponds 
to Q ^ at all temperatures < t < 1. In this case, 
the phase transition occurs into a single plane wave, or 
helicalfiS superconducting state with rj{r) = ryoe*''^. In 
particular, at weak fields near i = 1, the maximum of 
I(y, z) is achieved at yc = ~bz, which corresponds to 
Q = ~5h. For the critical field one obtains: 

At low temperatures, we use the fact that maxj, X(y, z) = 
— (p+/2) ln(8e'''~-'^z) in the limit z — > 00 (for < /9_). 
Therefore, 

~^\^-^ = ^.^''-'''^- (55) 

In these expressions we have omitted the terms propor- 
tional to . The temperature dependence of the critical 
field is shown in Fig. [51 both in the LOFF state, for 
p+ = /3_ = 1, and in the helical state, for = 0.8, 
p_ = 1.2. 

The most notable feature of the phase diagram is a 
considerable weakening of the paramagnetic pair break- 
ing at low temperatures, which is manifested in the diver- 
gence of the critical field at f — > 0."^^ Similar behavior has 
also been found in two-dimensional noncentrosymmetric 
superconductors li^ii^ 

We would like to mention that the paramagnetic pair 
breaking disappears altogether in the (rather unrealistic) 
extreme single-band case \S\ = 1, which corresponds to 
a very strong SO coupling, g ^ 00, see Eq. ([^. Al- 
though a quantitative treatment of this case is beyond 
the limits of applicability of our model, one can use sim- 
ple arguments to show that the critical temperature is 
not affected by the applied field. We first note that the 
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FIG. 2: Paramagnetic critical field h = ah in the clean case, 
for the LOFF state {S = 0) and the helical state with S = 
—0.2. Inset shows the temperature dependence of the wave 
vector Q of the order parameter modulation along the field. 

Zeeman interaction causes an anisotropic deformation 
of the electron bands: £,\{k) + \{g /2)pBkH . 

The intraband pairing of electrons with opposite mo- 
menta costs additional energy at TJ 7^ 0, resulting in 
the uniform superconductivity being suppressed by the 
Zeeman field. On the other hand, it is easy to see that 
the field-induced deformation of the bands amounts to 
shifting the bands in the opposite directions along the 
field: ^\{k) £_x{k + q\), where qx = \{g /2)p,bH /vp- 
If the Cooper pairs in the and "— " bands are com- 
pletely decoupled, or if there is just one band present, 
then the band shifts can be eliminated by independent 
gauge transformations, and the Zeeman pair breaking 
will be absent. According to Ref. '23', in the BCS-contact 
model ([ni the pairs in the two bands are in fact strongly 
coupled, so that both condensates are characterized by 
the same wave function rj{r). The band shifts cannot be 
eliminated simultaneously in both bands by any gauge 
transformation. Therefore, in general there is some en- 
ergy penalty associated with the Cooper pairing (both 
uniform and nonuniform) at ^ 0, compared with the 
zero-field case. 



B. Disordered case 

At arbitrary disorder, the paramagnetic critical field 
is obtained by numerical solution of Eq. (|38p . with Wn 
given by expressions ([5T|l . At (5 = 0, we find that the 
LOFF modulation is suppressed by impurities, see Fig. 
[21 and disappears at C > — 1-16. In contrast, the he- 
lical modulation at (5 ^ survives even if the impurity 
scattering is strong, as shown in Fig. [H In both cases, 
the low-temperature divergence of the critical field char- 
acteristic of a clean system, see Eq. ([55]) , is removed by 
disorder. 

In the "dirty" limit C S> 1 , one can obtain an equation 
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FIG. 3: Paramagnetic critical field in the LOFF state {5 — 
0) for different strengths of disorder. Inset shows the order 



The pair-breaker strength is characterized hy a = /QC,. 
Analytical expressions for the critical field can be ob- 
tained in the weak-field limit: 

h\t^i=\m{i-tf'\ (57) 

V 

and also at low temperatures: 

^1-0 = (58) 

Thus we come to the conclusion that, while the critical 
field near Tco is enhanced by disorder, the impurity re- 
sponse at low temperatures is non-monotonic: Initially, 
the disorder cuts off the singularity at T ^ 0, thus reduc- 
ing the critical field. However, as the disorder increases 
one eventually reaches the universal regime, in which the 
critical field grows as -/C,. 




FIG. 4: Paramagnetic critical field in the helical state with 
5 = —0.2, for different strengths of disorder. Inset shows the 
order parameter modulation along the field. 



for the critical field in a closed form. We shall see that 
both Q and h scale as ^ C,. Therefore, one can 
perform the Taylor expansion in Eq. (j51[) : 



Wn[Q) 



1 



{2n+l)t + C 



1 



W 



{{2n+l)t + Cf 



where W — {Q^ + /i^ + 2SQh)/3. The main contribution 
to the Matsubara sum in Eq. ([55]) comes from {2n+l)t ^ 
C, and we obtain: ln(l/t) = 4'(l/2 + W/2Ct) - *(l/2). 
It is easy to see that the maximum critical temperature 
is achieved when W has a minimum with respect to Q, 
which happens for Q — —Sh, at all temperatures. Thus 
we arrive at a well-known universal equation, which de- 
scribes magnetic pair breaking in superconductors)^ 



lni = *fi + ^V*f- 
t \2 t \2 



(56) 



C. First-order phase transition 

Let us now compare the paramagnetic critical fields 
found above with the critical field of a first-order phase 
transition (FOPT) into a uniform superconducting state 
at zero temperature. For simplicity, we neglect the dif- 
ference between the band densities of states and set 
g — 2. Following the arguments of Clogston and 
Chandrasekharj— we obtain: ^ibHfopt = ('7o/'\/2)(l — 
Xs/xp)"^^^, where r\Q = (7r/e'^)Tco is the gap magni- 
tude at zero temperature and field, Xs is the residual 
susceptibility at T = in the superconducting state, and 
Xp = 2jj,^Np is the Pauli susceptibility in the normal 
state. Therefore, 



1 



V2f 



XpJ 



-1/2 



(59) 



According to Ref. the residual susceptibility in our 
model is given by the following expression: 



XP 3 3 



(60) 



where x — 2r/377o, and 



$(x) 



TT 

2x 



7r\/r 



: arctan 



[at a; > 1 this function is evaluated using arctan(ia;) = 
itanh~^(a;)]. In the clean limit, <&(0) = and Xs/xp = 
2/3. Therefore, hpoPT = ^/ij2e-^ ~ 0.69, which is 
much lower than the divergent expression (j55p . In the 
"dirty" limit x oo, we have $(x) ~ 1 — 7r/2a;, and it fol- 
lows from Eq. ^ that Jifopt = ^/2CJt^. Although 
the critical field of the first-order phase transition into a 
uniform superconducting state increases with disorder, it 
remains lower than the critical field of the second-order 
phase transition into the helical state, which is given by 
the expression ([51 



GENERAL CASE: ORBITAL EFFECTS 



Solution of Eq. ((38)) at all temperatures and for ar- 
bitrary values of a and ^ can only be obtained numeri- 
cally. To facilitate numerical analysis, we represent the 
integrals in Eq. (|39p in a somewhat shorter form by in- 
troducing the Cartesian coordinates p = (pi , p2 , Ps ) j such 
that p = IpI, ps = ps, and — s^) — p\ 
Using the Fourier transform 



drpe P — — = — arctan- 



pI 



p^ k 



we obtain 



1 



k dk 



X arctan 



^fc2 + (Q + XahY 

lN{k), (61) 



(2n-Kl)i + C 



where 



= ^/ dVxe-'^^-e-''''i/4i^.(^). (62) 

In the purely paramagnetic limit, /jy = 27r(5(fc), and one 
recovers Eq. (I5T|) . 

The maximum critical field corresponds to TV = 0, in 
which case /o(fc) = (2//i) exp(— /c^/ft,). The order pa- 
rameter is nonuniform along the direction of the applied 
field, while its dependence on the transverse coordinates 
is given by the usual Abrikosov vortex solution. As in 
the purely paramagnetic limit, the phase diagram turns 
out to be different for 5 = (LOFF vortex state) and 
(5 7^ (helical vortex state). 

Solution of Eqs. ^ and ^ shows that the LOFF 
state is suppressed by the orbital effects, even in the 
clean case. There is a critical strength of the param- 
agnetic interaction, ac, below which the LOFF state dis- 
appears, i.e. Q = at all temperatures (in the clean limit 
Q!c — 0.66). This is qualitatively similar to the way the 
orbital interaction affects the LOFF state in centrosym- 
metric superconductors, see Ref. [s^ . In contrast, the 
helical modulation at (5 7^ 0, albeit weakened by the or- 
bital effects, does not completely disappear until a = 0. 

In order to study the combined effect of the orbital 
interaction and impurities, we focus on the helical vortex 
state with 5 — —0.2. We consider a = 0.2 [which is 
close to the estimates for the Li2(Pdi_a:,Pta;)3B family 
of superconductors, as obtained from Eq. (j43p using the 
data from Ref. and also a = 2.0. The upper critical 
field curves for a range of disorder strengths are shown 
in Figs. [5] and [SI We see that if a is not too large, 
disorder produces a monotonic enhancement of Hc2 at 
all temperatures. 

In the "dirty" limit, it is again possible to obtain a rel- 
atively simple equation for the upper critical field. Re- 
peating the arguments from Sec. IIVBI one can show 
that that equation has the universal form (|56p . with the 
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FIG. 6: Upper critical field hc2{t) in the helical vortex state 
with a = 2.0 and 5 = —0.2, for different strengths of disorder. 



paramagnetic and orbital interactions both contributing 
to the pair-breaking strength: 



2u2 



h + ah 
~~6C 



(63) 



The maximum critical field is achieved for Q = —6ah. 
Near Tco, superconductivity is suppressed by magnetic 
field according to Eq. (|^ . in which ai = ■k^/12( + 
OiC^) and 02 = (7r2/12C)a2 + 0{C^). At low temper- 
atures, returning to dimensional variables, we arrive at 
the following expression for the upper critical field: 




(64) 



In the weak SO coupling limit, corrections to this expres- 
sion are of the order of 5'^. In the orbital limit a ^ 0, one 
obtains Hc2{0) = {3/2ne'^){r /Tco), while in the param- 
agnetic limit a 00, Eq. ([55)1 is recovered. We see that 
the upper critical field is enhanced by disorder. 
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FIG. 7: Ratio Q/{~Sah) in the helical vortex state with S — 
—0.2, for the following disorder strengths: ^ = 0.0, 0.5, and 
5.0 (from top to bottom). 

As for the wave vector of the helical modulation, the 
expression (|45p is exact only in the Ginzburg-Landau 
regime and in the "dirty" limit. In clean systems at 
low temperatures, deviations of Q / (—Sah) from unity 
become quite substantial, see Fig. [71 

VI. CONCLUSIONS 

We have calculated the upper critical field in non- 
centrosymmetric superconductors at arbitrary temper- 
atures, both in the clean case and in the presence of 
scalar impurities, using as an example a crystal of cubic 
symmetry with a BCS-contact pairing interaction. Our 
derivation relies on the assumption that the SO band 
splitting is small compared to the Fermi energy. Then 
one can neglect the disorder-induced triplet component of 
the gap function, which considerably simplifies the calcu- 
lations, without losing information about the interesting 
physics of the problem, in particular about the properties 
of nonuniform superconducting states. 

As in the centrosymmetric case, the pair breaking is 
provided by both the orbital and the Zeeman (paramag- 
netic) interactions. If the former is dominant, i.e. a — *■ 0, 
see Eqs. (pij) and then the lack of inversion sym- 

metry does not have any appreciable effect on the upper 
critical field, which is described by the standard Helfand- 
Werthamer theory. However, as a increases, so do the de- 
viations from the centrosymmetric case. In the extreme 
paramagnetic limit, corresponding to a ^ oo, the critical 
field diverges at T in the clean case. 

Impurities suppress both the orbital and paramagnetic 
pair breaking, see Eqs. ([5^)1 and ([SH]) . resulting in an en- 
hancement of the upper critical field. In the strongly 
paramagnetic limit, the effects of disorder on the low- 
temperature critical field are nonmonotonic: At weak dis- 
order Hc2 decreases, but in the "dirty" limit the overall 
enhancement of the critical field takes over. 



The spatial structure of the superconducting order pa- 
rameter is, in general, different from the centrosymmetric 
case. As soon as the difference between the densities of 
states in the SO-split bands, and N^, is taken into 
account, the system exhibits a helical instability, in which 
the order parameter has a phase gradient along the ap- 
plied magnetic field (in a cubic crystal). The helical su- 
perconducting state turns out to be robust with respect 
to both the orbital pair breaking and disorder. In par- 
ticular, the wave vector of the superconducting instabil- 
ity in in the Ginzburg-Landau regime and in the "dirty" 
limit is given by Qz — —26a{H/Ho)^Q^. In contrast, the 
LOFF state, which occupies the low-temperature part of 
the phase diagram at = , is quickly suppressed by 
both the orbital interaction and disorder, similarly to its 
counterpart in the centrosymmetric case.'^*''^^ 

Let us discuss the application of our results to 
the family of cubic noncentrosymmetric compounds 
Li2(Pdi_a;,Pta;)3B, where x ranges from to 1. The crit- 
ical temperature varies from 7-8 K for a; = to 2.2-2.8 
K for X = 1. The electronic band structure also ex- 
hibits considerable variation: The SO band splitting is 
strongly anisotropic, and can be as large as 30 meV in 
LiaPdgB, and 200 meV in LiaPtgB (Ref. E^. Experi- 
mental data^i^^ seem to agree that Li2Pd3B is a con- 
ventional BCS-like superconductor with no gap nodes. 
In contrast, the gap structure in Li2Pt3B is still a sub- 
ject of controversy. While earlier experiments, see Refs. 
ISTI . suggested the presence of lines of nodes in the gap, 
the recent /xSR and specific heat data^S, have found no 
evidence of those. Moreover, according to Ref. [H, the 
whole Li2(Pdi_a;,Pta;)3B family of compounds are single- 
gap isotropic superconductors. If this is the case then 
our model based on a BCS-contact pairing Hamiltonian 
should be applicable. In these compounds, the param- 
agnetic effects seem to be rather weak: the parameter 
a varies from 0.10 for x = to 0.15 for x = 1 (us- 
ing the data from Ref. [ssh . This explains why the ex- 
perimental Hc2 curves in these materials show a good 
agreement with the Helfand-Werthamer theory. Still, 
the absence of inversion symmetry should manifest it- 
self in a long-wavelength helical modulation of the or- 
der parameter along the applied field. Using the maxi- 
mum values of the SO band splitting from Ref. [36l . we 
have: Eso - 30 meV in Li2Pd3B, and Eso - 200 meV 
in Li2Pt3B. Then one can obtain the following order- 
of-magnitude estimates for the wave vector of the heli- 
cal modulation: Qz^q ~ 10~^H/Hc2{0) in Li2Pd3B, and 

Qzio - 10-2i?/iJ,2(0) in Li2Pt3B. 
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APPENDIX A: SPECTRUM OF ^00(0^™) 

The problem of finding the upper critical field is re- 
duced to the calculation of the average normal-state 
Green's function in a uniform magnetic field. The differ- 
ence from the centrosymmetric case is due to the presence 
of the SO coupling term in the Hamiltonian ([5]). Follow- 
ing Ref. 139|, we represent the Green's function before 
disorder averaging in a factorized form: 

Go.p{ri,r2;iOn) = Gc.0{ri,r2;uj^)e''^(-^''^'\ (Al) 

where ip{ri,r2) — (e/c) Jl^^ A{r)dr, and the integration 
is performed along a straight line connecting ri and r2. 
From Eq. p2|) we obtain an equation for the gauge- 
invariant factor: 



(iw„ - hi)G{ri,r2;LJn) ^ S{ri - r2), (A2) 



where 



hi = --^+-foKia+^l^iBH& + Uiri)-eF, (A3) 
Zm* Z 

and Ki = -iVi + {e/ 2c)H x (n - r2). 

As discussed in Sec. Ill Al in the limit of weak SO cou- 
pling it is sufficient to consider only the singlet-singlet 
terms in the Cooper impurity ladder. After disorder av- 
eraging, the kernel ([5^ takes the following form: 



3^00(^1, 'r2;a;„) = 3^00(^1 - r2, w„)e 



2iip{ri ,r2) 



(A4) 



where the translationally-invariant part is defined by its 
Fourier transform as follows: 



3^00 



1 



2ttNf 



(2^ 



tr [G{k + q, 0Jn)^2G'^ {-k, -ujn)^2]- (A5) 



The Green's functions here are disorder-averaged solu- 
tions of Eq. (|A2p . The next step is to use the iden- 
tity e2^'^(''i'''2)77(r2) = e-^(''i-^2)^i?7(ri), where D = 
— iV + {2e/c)A, from which it follows that the opera- 
tor 3^00(^71) can be obtained from 3^oo(9: ^n) by replacing 
q ^ D. The gauge-invariant part of the Green's function 
contains the effects of the SO coupling, the Zeeman in- 
teraction, the orbital Landau quantization, and disorder. 
We treat the magnetic field effects on G perturbatively, 
which is legitimate if the Zeeman energy is small com- 
pared with the SO band splitting, i.e. {g/2)fiBH ^ Eso, 
and also if the temperature is not very low, so that the 
Landau level quantization can be neglected. 

In the clean case, keeping only the terms linear in H 
in Eq. I|A2|) . we have 



fiQ + Hm)G{r,L 



Sir), 



(A6) 



where ho = eo(— iV) +7(— iV)o- is the zero-field Hamil- 
tonian, and rh = —^B[{g/'2)o- + (R x tt)], with tt = 



{m/m*){—i'V) -|-m7o6", has the meaning of the magnetic 
moment operator for band electrons. The eigenvalues of 
ho are given by ^\{k), see Eqs. It is straightfor- 

ward to show that the linear in H effects on G are due 
to the Zeeman interaction (the first term in rh). Then, 

G{k,uJn) can be represented in the form with the 

band Green's functions now given by 



Gx{k,Un) 



1 



Uk)-M.g/2)fiBkH 



(A7) 



Thus, the electron bands are deformed by the magnetic 
field: Uk) ^ Uk) + X{g/2)tiBkH. 

In the presence of impurities, again keeping the mag- 
netic field only in the Zeeman term in Eq. (|A3p . we 
obtain for the average Green's function: 

G{k,Un) = [icj„ - eo(fc) - 7(fc)o- - Simp(t^„)]"\ (A8) 
where 7(fe) = 70^ + {g/2)fiBH, and 

d^k 



G{k,uj^) (A9) 



is the impurity self-energy in the self-consistent Born 
approximation. We seek the self-energy matrix in a 
spin-diagonal form: Yiimp{ijJn) — —'iPi^n)o'Q- Then, 
G{k,ujn) = J2\ ti\ik)G\{k,ujn), where a)„ = w„-|-/?(a;„), 
and the band projection operators and the band Green's 
function are given by expressions (fH|) and (|A7|) respec- 
tively. Substituting this into Eq. (jA9p and neglecting the 
energy dependence of the band densities of states on the 
scale of the Zeeman energy, we find: cun = Ci^„ -f F sign . 
Therefore, 



(5(fc,c^„) = ^nA(fc) 



1 



(AlO) 



i^^n - ^A(fc) - Hg/2)^iBkH + iT signcjn 
Inserting the last expression in Eq. (|A5|1 . we obtain: 



PA 



1 



F + zJIa signa;„ 



(All) 



where p\ = N\/Nf, and nx{k,q) ~ VFk{q + 
Xg^sH / vf) /2. Next we use in Eq. (|A11[) the identity 
a^^ — due~°'^'' and make the substitution q ^ D 

to represent 3^00 ("-^n) as a differential operator of infinite 
order: 



3^00 (w„) 



Here 



Ox - ( exp 



due-"(l""l+r'^PAOA. (A12) 



*™^.fcfr> + A^) ) (A13) 
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do not depend on a;„, because the zero- field Fermi sur- 
faces are invariant under k —k. 

To find the eigenfunctions and eigenvalues of Ox, we 
follow the procedure outlined in Ref. We choose the 
z-axis along the external field, so that H = Hz, and 
introduce the operators 

^^^^^D^^iDy^ a3=^/^i5., (A14) 

where £h = \/c/eH is the magnetic length. It is easy to 
check that a+ = a[_ and [a_,a+] = 1, therefore a± have 
the meaning of the raising and lowering operators, while 
03 — flg commutes with both of them: [03, a±] = 0. We 
use the basis of states \N,p) (Landau levels), such that 

a+\N,p) = VN+l\N + l,p) 
a^\N,p) = ^/N\N-1,p) 
a3\N,p) =p\N,p), 



where iV = 0, 1, and p is a real number. Then, 

dx\N,p) = - [ d9 sine e-'^'^P+^Po^"'"''^ 

Jo 27r 
2 J-i 

xLn[v^(1-s^)]\N,p), (A15) 

where v = {vp /'^■Ih)u, po = gfiBHln /vp, and Ln{x) 
is the Laguerre polynomial of degree N . We see that 
the Landau levels \N,p) are eigenfunctions of 0\ and, 
therefore, of 3^oo('^n)- Introducing 5 = {p+ — P-)/2 and 
summing over A in Eq. (jA12p , we arrive at the expression 
0. 
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